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The majority of dynamical studies in power systems focus on the high voltage transmission grids where 
models consider large generators interacting with crude aggregations of individual small loads. However, new 
phenomena have been observed indicating that the spatial distribution of collective, nonlinear contribution of 
these small loads in the low- voltage distribution grid is crucial to outcome of these dynamical transients. To elu- 
cidate the phenomenon, we study the dynamics of voltage and power flows in a spatially-extended distribution 
feeder (circuit) connecting many asynchronous induction motors and discover that this relatively simple 1+1 
(space+time) dimensional system exhibits a plethora of nontrivial spatio-temporal effects, some of which may 
be dangerous for power system stability. Long-range motor-motor interactions mediated by circuit voltage and 
electrical power flows result in coexistence and segregation of spatially-extended phases defined by individual 
motor states-a "normal" state where the motors' mechanical (rotation) frequency is slightly smaller than the 
nominal frequency of the basic AC flows and a "stalled" state where the mechanical frequency is small. Transi- 
tions between the two states can be initiated by a perturbation of the voltage or base frequency at the head of the 
distribution feeder. Such behavior is typical of first-order phase transitions in physics, and this 1+1 dimensional 
model shows many other properties of a first-order phase transition with the spatial distribution of the motors' 
mechanical frequency playing the role of the order parameter. In particular we observe (a) propagation of the 
phase-transition front with the constant speed (in very long feeders); and (b) hysteresis in transitions between 
the normal and stalled (or partially stalled) phases. 



Popular Summary: Large electrical generators interacting 
over national- scale electrical transmission grids constitute a 
well- studied dynamical system. Rarely discussed and poorly 
understood are the dynamics of neighborhood-scale distribu- 
tion grids extending from transmission substations to the mul- 
titude of individual customers. However, the changing na- 
ture of electrical loads, e.g. the increasing prevalence of in- 
duction motors in residential air conditioning units, is creat- 
ing distribution-grid dynamical processes that, when excited 
by unremarkable transmission-grid disturbances, lead to ir- 
reversible transitions with major impact on the reliability of 
transmission grids. Here, we present a unique model and anal- 
ysis of these dynamics that allow analogy with other physical 
systems enabling rapid progress by leveraging knowledge de- 
veloped in physics. 

In contrast to usual approaches, we develop a spatially- 
continuous model of distribution-grid dynamics to investi- 
gate the collective dynamics that arise when many induc- 
tion motors, which are individually nonlinear and hysteretic 
(bi-stable), are coupled via power flows and voltage evolu- 
tion within a distribution grid. Normal- size perturbations ex- 
cite these collective dynamics initiating soliton-like fronts that 
travel though the distribution grids where passage of the fronts 
results in transitions of individual motors between a normal 
state and an undesirable stalled state. Individual bi- stability of 
each motor promotes globally hysteric behavior that is remi- 
niscent of well-known first-order phase transition dynamics 
found in other physical systems. 

Important extensions of this physics-based understanding 
include the ability to model the dynamics of the billions ac- 
tive "smart grid" loads predicted to revolutionize the electrical 
power system. 



I. INTRODUCTION 

Power systems are used to generate and transfer energy 
to electrical loads. In today's grid, generation is primarily 
done at large, centralized power stations (~100's of MW) 
and the transfer primarily occurs via alternating currents (AC) 
in national-scale, highly-meshed, high-voltage transmission 
grids. A subset of nodes (substations) in the transmission 
grid transform the high voltage to a medium voltage level and 
interface to distribution grids, however, another change oc- 
curs at the substation. The meshed network of national- scale 
transmission changes to many radial or tree-like structures in 
the distribution system whose spatial extent is only ~ 1-10 
km. Each radial circuit, also called a "feeder", distributes the 
power delivered to the substation by the transmission system 
to the thousands of small electrical loads (~ 1 kW) spatially 
spread along its length. 

Even though AC electrical generation and transmission 
grids are extended over large spatial scales, they are syn- 
chronized, i.e. power flows over the transmission lines cre- 
ate dynamical coupling between the large rotating generators 
forcing them to rotating in unison. Perturbations to this syn- 
chronized system results in dynamics and transients spanning 
a large range of temporal scales from milliseconds to many 
minutes. Many studies have addressed the dynamics of these 
large-scale transmission grids. Although many unresolved dy- 
namic problems in the transmission grid remain, recent years 
have witnessed new phenomena that have refocused our atten- 
tion on dynamics and transients occurring in the smaller-scale 
distribution circuits 1 2 . Although these phenomena occur on 
smaller scales, they involve collective behavior of many indi- 
vidual small nonlinear loads, and even coupling several dis- 



2 



tribution circuits, creating a significant impact on the larger- 
scale transmission system. 

Perhaps the most drastic of these phenomena is voltage 
collapse^. Here, a quasi-static increase in loading pushes the 
distribution feeder to a bifurcation where the stationary nor- 
mal/high voltage solution is lost and the feeder "collapses" 
to an undesirable low voltage solution that is dangerous for 
power system stability. However, even for moderately loaded 
feeders that are far from this critical point, the nonlinearity of 
electrical loads may result in the emergence of multiple sta- 
tionary solutions. In contrast to the case just described, these 
solutions cannot be reached via quasi- static evolution of the 
electrical loads. Instead, a significant and nonlinear pertur- 
bation to the feeder creates a dynamical trajectory that ter- 
minates in one of these additional solutions that may also be 
a "bad solution" from the standpoint of voltage level, power 
system losses, stability, and equipment damage. To under- 
stand the possibility of this unwelcome outcome one needs to 
go beyond the traditional static description and analyze dy- 
namics of the distribution system. 

The spatio-temporal dynamics we seek to describe occur 
within an individual radial distribution feeder that connects 
many (~ thousands) small loads to a substation. We con- 
sider electro-mechanical dynamics occurring on scales rang- 
ing from fractions of a second to tens of seconds and analyze 
the spatial distribution of power flows along the circuit and 
the spatio-temporal transients stimulated by exogenous distur- 
bances in the voltage and base frequency at the head of the dis- 
tribution feeder. Such disturbances, which primarily originate 
from faults and/or irregularities in the high- voltage transmis- 
sion system, propagate through from load to load via power 
flows in the distribution feeder. The propagation is affected 
by the electro -mechanical response of individual loads, and 
here we focus on the effects of nonlinear loads such as asyn- 
chronous (i.e. induction) motors. Typical electro -dynamic 
transients propagate with speed comparable to the speed of 
light and damp out in tens of milliseconds and thus are not 
important in our analysis. On the other hand, composition 
of loads connected to a distribution feeder changes on much 
longer time scales (~ minutes) and are taken as fixed on the 
time scale of electro-mechanical dynamics of interest in this 
work. 

An interesting and key feature of the dynamics is the long- 
range coupling between the spatially distributed loads created 
by power flows along the feeder. These coupled dynamics are 
nontrivial to model and investigate, however, they are also ex- 
tremely important for practical power engineering because our 
model solutions reveal serious problems for control and oper- 
ation of power systems. One such problem that motivates our 
study is the phenomeno n of the Fault-Induced Delayed Volt- 
age Recovery (FIDVR) A FIDVR event is typically 
initiated by a fault on the transmission grid near a substation 
creating large fault currents that temporarily depress the volt- 
age at the substation, perhaps for as little at two cycles of the 
nominally 50/60 Hz AC frequency (~ 30 msec). The voltage 
depression propagates into the substation's distribution feed- 
ers causing an almost instantaneous reduction in the electrical 
torque generated by the connected induction motors, however, 



the mechanical torque on the motors does not change instan- 
taneously and the motors begin to decelerate. If the transmis- 
sion fault and voltage depression last long enough, many of 
the induction motors along the feeder may stall. When the 
fault on the transmission grid is cleared, the voltage at the 
substation returns to near normal levels, however, a stalled in- 
duction motor draws large in-rush currents while at near zero 
rotation speed (mechanical frequency). The time synchroniza- 
tion of these in-rush currents cause large voltage drops in the 
distribution feeder and may hold the voltage at locations re- 
mote from the substation below a critical voltage for restart- 
ing. Crucially, these remote motors remain stalled (near zero 
mechanical frequency), and their large current draw stabilizes 
this spatially-extended, partially stalled state. 

The tendency for a transmission fault to result in FIDVR 
may depend on many fault, distribution feeder, and motor pa- 
rameters, e.g: voltage drop magnitude and duration; length, 
resistance, and reactance of the feeder; type, rotational inertia, 
and density of the induction motor loading; and possible post- 
fault corrective control actions. The qualitative description 
of FIDVR given above provides an intuitive understanding 
of how some of these parameters affect the dynamics leading 
to the undesirable and potentially dangerous partially- stalled 
state. It certainly indicates that, without some sort of cor- 
rective actions, FIDVR will become more and more frequent 
because of the recent trends to more air conditioning driven 
by easy-to-stall, low inertia motors. However, the dynam- 
ics that lead to FIDVR are not generally understood and thus 
presumed somewhat mysterious in power engineering prac- 
tice. The goal of this manuscript is to provide understanding 
of these interesting and practically important distribution grid 
dynamics. 

FIDVR is an example of a broader class of problems where 
physics and dynamical systems modeling can provide signif- 
icant insight and predictive power that are generally lacking. 
Another example is given by electro-mechanical waves propa- 
gating through transmission grids 12 13 , as well as tra nsien ts as- 
sociated with the loss of synchrony in power systems-SEI] The 
key unifying feature of all these phenomena is in the nontrivial 
interplay of spatial coupling of individual (possibly nonlinear) 
dynamics via power flows over the electrical network. We be- 
lieve that important insights into these complex dynamics can 
be gained by approaching such spatio-temporal phenomena 
from a homogenized prospective, i.e. studying the electrical 
grid not as a set of individual devices but rather as a spatially 
extended and continuous medium in the limit where the num- 
ber of individual elements of the power system becomes in- 
finite. This abstract continuous -medium approach, pioneered 
for the case of electro-mechanical waves over transmission 
systems irPEU anc [ f or th e case f a radial/linear distribution 
system irP^is advantageous as it enables (a) a simpler analysis 
and deep qualitative physical understanding of the underlying 
phenomena (e.g. of FIDVR and electro-mechanical waves); 
(b) flexibility in simulations; and (c) developing model reduc- 
tion algorithms for faster state estimation and system simula- 
tion. 

Motivated by the discussion above, the main goal of this 
manuscript is to formulate the simplest but still realistic 1+1 
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(space+time continuous) model that predicts and explains in- 
teresting and important spatio-temporal phenomena in a dis- 
tribution feeder loaded with induction motors which can be in 
a normal or stalled state. The most important results reported 
in this manuscript are 



• Extending the previous works^EIHIS we show that 
if the local voltage falls sufficiently low, an individ- 
ual asynchronous motor can be in either of the fol- 
lowing two states: (1) a normal state characterized by 
mechanical frequency uj which is slightly lower than 
the base electrical frequency ujq of the system; (2) a 
stalled state characterized by low or zero mechanical 
frequency. Both states are locally stable but can evolve 
into each other under sufficiently large perturbations. 

• In sufficiently long feeders, a partially stalled phase can 
emerge where the feeder splits into head and tail parts 
with the motors of the head (tail) being in the normal 
(stalled) state. Motors in the stalled state may occupy 
the entire distribution feeder or, if the feeder is very 
long, stalled portion can co-exist with the normally run- 
ning one. In the latter case, there exist multiple, par- 
tially stalled phases characterized by different propor- 
tions of the head (normal) and tail (stalled) parts. 

• The steady partially stalled phases can be interpreted as 
showing coexistence of the two states where the local 
mechanical frequency (of the motors) play the role of 
"order parameter". Transitions between the phases are 
classified as first order, using standard physics termi- 
nology. 

• These transitions are hysteretic (not reversible), i.e. a 
perturbation leading to transition from the normal phase 
to the stalled phase is not an inverse of the perturba- 
tion leading from the stalled phase to the normal phase. 
The dynamics of the two transitions are also different, in 
particular fronts of the phase transitions have different 
shapes, and to stabilize one transition can take signifi- 
cantly longer than the other. 

Material in the manuscript is organized as follows. Static 
models of single induction motor, two-bus system and the 
DistFlow equations discussed in Sections |ll A|lIB 

and III A respectively. A dynamic model of a distribution 
feeder loaded with induction motors and 1+1 space-time con- 
tinuous model of this feeder are introduced and discussed in 



Sections IIIB and IIIC Section IV A discusses how single- 
motor bi- stability translates into the emergence of multiple 
phases of the feeder with supporting numerical experiments 
in Section |IVB] and special features of the phase transitions 
Section |V| is devoted to in-depth discussion 
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in Section 

of the results of our numerical experiments. Section |VA| dis- 
cusses the dynamics following a fault at the head of the feeder. 
Sec tion |VB| analyzes the recovery from a stalled state. Sec- 
tion [VCfexplores the phase space of parameters that governs 



whether or not a feeder will enter a stalled or normal state fol- 
lowing a fault. Finally, we summarize and describe a path for- 
ward in Section [Vl] Auxiliary information explaining details 



of our simulations can be found in Appendix [A] Appendix 
[B] provides captions for the illustrative movies of the phase 
transition and fault recovery processes available as a Supple- 
mentary Information (SI) to the manuscript. 



II. SINGLE INDUCTION MOTOR MODELS 

A. Static Motor Model 

Induction motors play a significant role in FIDVR, as it is 
currently understood. Here, we describe the features of in- 
ductions motors that are important for the rest of our work. 
Although motor dynamics will be added later, we first adopt a 
simple static electrical model of an induction (asynchronous) 
motor rotating at mechanical frequency uj and connected to a 
distribution circuit being driven at a base frequency ujq p2 
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Here, P and Q are real and reactive powers drawn by the 
motor; s = 1 — uj/ujq is the slip parameter of the motor, 
< s < 1; v is the voltage at the motor terminals; X m , R m 
are internal reactance and resistance of the motor, usually 
Rm/Xjn = 0.1 -T- 0.5. For steady rotation frequency at uj, 
the balance of electric and mechanical torques for the induc- 
tion motor is 



UJ 



= T(uj/uj ), 



(3) 



where T(uj/ujq) is the rotation speed-dependent torque ap- 
plied to the motor shaft by the mechanical load, which is typ- 
ically parameterized by 



T(uj/uj ) = T 



UJ 



(4) 



Here, T is a reference mechanical torque and a is indica- 
tive of different types of mechanical loads with a = 1 typical 
of fan loads and a < 1 typical of air-conditioning loads. If 
a < a c ~ 1 and To is fixed, one observes the emergence 
of three solutions when v is in a range between two spinodal 
voltages v~ and v+, i.e. for v~ < v < v+ (see Fig. [I]). These 
solutions have widely different uj which leads to hysteresis 
and the interesting dynamical behavior explored in the rest of 
this manuscript. 



Stability analysis (see Section II B for details) shows that 
the two extreme solutions (uj ~ and uj w ujq) are both stable 
while the solution in the middle is unstable. The consequence 
is hysteretic behavior of the motor frequency uj as a function 
of the voltage v, as displayed in Fig. [2] Starting in the high- 
voltage normal state (say v ~ 1), we decrease v slowly along 
the dashed red curve passing through state d. If we further 
decrease v to state c, the normal state suddenly disappears, and 
the motor makes a transition to the stalled state at a. Similarly, 
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FIG. 1 : Electric and mechanical torques as functions of the mechan- 
ical frequency cj/cjo for a range of motor terminal voltages v, refer- 
ence mechanical torque To = 0.32, and a = 0.1. For the v = 0.86 
electrical torque curve (red), there are three equilibrium solutions in- 
dicated by intersections with the mechanical torque curve (black). 
The solution with the highest u/uo is the "normal" stable solution 
with the induction motor rotating near the grid frequency cjo- The 
"stalled- state" with uj/ujq ~ is also stable while the intermedi- 
ate solution is unstable. For v > = 0.9 (light blue curve), 
there is only one solution corresponding to the normal state. For 
v < v~ = 0.83 (green curve), there is only one solution correspond- 
ing to the stalled state. The points (a, b, c, d) correspond to the same 
labels in Fig. [2] 



if we start from the low- voltage stalled state (say v ~ 0.75 on 
the black curve) and v is increased slowly through state a to 
b, the stalled state disappears and the motor makes a transition 
to the normal state at d. The states (a, 6, c, d) are also marked 
in Fig.[T] and the same hysteresis loop can be traced out there. 

After elimination of the auxiliary variable u, Eqs. ( |1|2|3|4| ) 
describe the dependence of the power flows (P,Q) on 
terminal- voltage. If the reactance-to-resistance ratio of the 
motor, Xm/Rm, is sufficiently small, the hysteresis observed 
for mechanical frequency in Fig. [2] translates into hysteresis of 
real and reactive powers as seen by the multi- valued depen- 
dence of (P, Q) on v in Fig. [3] Between the spinodal voltages 
v ~ and v+ , there exists three solutions with different values of 
(P, Q) for the same value of voltage. Of the three solutions, 
the top and bottom are stable while the middle solution is un- 
stable. Similar to Fig.[2j we can follow an adiabatic evolution 
of the motor terminal voltage. Following the reactive power 
curve (red) and starting in the normal state with v ~ 1, we 
decrease the voltage through state d and to state c. Any fur- 
ther reduction of v forces the motor to make a discontinuous 
jump to state a which is accompanied by a large increase in 
reactive power Q. Alternatively, we may start in the stalled 
state with v ~ 0.7 and slowly increase v through state a to 
b where the motor is forced to jump to state d accompanied 
by a discontinuous decrease in Q. As discussed later, these 
discontinuous jumps in Q play a significant role is stabilizing 
the spatially extended stalled state and in the recovery from 
the stalled state to normal state. 



FIG. 2: Hysteretic behavior of an induction motor in Fig. [T] as the 
voltage v at its terminals is varied. The dashed red (solid black) 
curves indicate the path of equilibrium states as the voltage v is 
decreased (increased) starting from the high-voltage normal (low- 
voltage stalled) state. The vertical lines at the spinodal voltages vf 
indicated the abrupt hysteretic transitions between states, vf corre- 
spond to the same labels in the legend of Fig. [T] and in Fig. [3] The 
points (a, 6, c, d) correspond to the states where the motor must make 
transitions from normal to stalled (c — >> a) and from stalled to normal 
(b^d). 
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FIG. 3: Typical (P, Q) versus v curves for the same motor as in 
Fig. [T] and [2] and described by Eqs. (1|2|3|4) . Three solutions (two 
stable and one unstable) are observed between the spinodal voltages, 
i.e. v~ < v < Vc . Solid and dashed curves (the latter partially cov- 
ered by solid) show trajectory of the system under adiabatic evolution 
starting from low-voltage and high- voltage regimes respectively. The 
points (a, 6, c, d) label the reactive power curve (red) and correspond 
to the same labels in Fig.[l]and Fig.[5] Here, X m /R m = 0.375. 



B. Dynamical Motor Model 

To study the important yet generic aspects of the dynam- 
ics in distribution circuits, we generalize Eqs. (T][3 ) to include 
induction motor dynamics. Considering for the moment a sin- 
gle induction motor, an imbalance in electrical and mechani- 
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cal torques will cause a change in the motor's rotational fre- 
quency given by 



*, d P 

M—uj — 

(If UJQ 



To 



Co) 



(5) 



where M is the motor's moment of inertia. Torque imbalances 
in Eq. §5§ can be driven in two ways: directly via changes in 
the base frequency ujq in Eq. p| or indirectly via changes P 
driven by changes in v in Eq. Q. The coupling of Eq. §5§ to 
ujo will not be strong because ujo is determined by the global 
balance of generation and load across the entire transmission 
system. We would not expect the local distribution dynam- 
ics under consideration here to affect ujo to a degree that we 
would have to consider its effect back on the distribution dy- 
namics via Eq. ([5]). Therefore, we can generally ignore the 
dynamics of uo and consider it an imposed exogenous param- 
eter. 

In contrast, changes in local voltage are strongly coupled to 
changes in the local flow of real and reactive power. As we 



have seen in Sec. |II A| changes in voltage can lead to dras- 
tic and hysteretic changes in a motor's frequency uo resulting 
in a strong coupling back to the dynamics in Eq. ([5]). There- 
fore, we must consider the possibility that voltage dynamics 
are important for distribution feeder dynamics. However, the 
dynamics of v are fundamentally different than uj because the 
relaxation of v is entirely electrical as opposed to the mechan- 
ical dynamics of u, i.e. 



Q 



s 2 X„ 



(6) 



where r is the characteristic time of this purely electrical pro- 
cess (See e.g. the Appendix of Pereira et apty. In our nu- 
merical experiments, we observe that the important basic phe- 
nomena discussed in the manuscript are much more sensitive 
to variations in the moment of inertia M than to variations in 
r and that setting r = still reveals the dynamical processes 
important for understanding FIDVR. With r = 0, Eq. ^ is 
now equivalent to its static version in Eq. §2§. 

With the dynamics now fully specified by Eq. ([5]), we can 
now justify the local (single motor) stability claims made in 
Section |n"A| by simple inspection of the equilibrium states in 
Fig. [T] Here stability/instability is understood in terms of the 
temporal decay/growth of small perturbations to the dynam- 
ics of Eq. ([5]). The black curve in this Fig. [T] represents the 
mechanical torque on the motor while the colored curves are 
the electrical torques (each representing a different v). Con- 
sider state d in Fig. [T] which is representative of the normal 
states with uj/uj ~ 1. If the motor speeds up slightly (moves 
to the right along the light blue curve), the mechanical torque 
becomes larger than the electrical torque and the motor de- 
celerates returning to d. If the motor slows slightly (moves 
left), the electrical torque becomes larger while the mechani- 
cal torque decreases returning the motor to state d. All of the 
normal states, i.e. those like d with uj/uq ~ 1, have similar 
behavior and are therefore stable. State a in Fig.[T]is represen- 
tative of the stalled states, and following the same logic, we 
find that the mechanical torque is larger for a higher uj (and 



smaller for a lower uj) showing that all of the stalled states 
are stable. Following the same logic, we find the the states in 
between the normal and stalled states (e.g. the state given the 
intersection of red and black curves in Fig.[T]) are unstable. 



III. SPATIALLY-CONTINUOUS FEEDER POWER FLOW 
MODEL 

In Section [TTJ we described the dynamics of isolated induc- 
tion motors, i.e. motors whose terminal voltage v is specified 
and not determined in part by interactions with other induction 
motors or electrical loads. In this Section, we consider power 
flow models responsible for creating the long-range coupling 
between the individual, local induction motor dynamics. We 
start with a well-known discrete power flow model which we 
then homogenize into a spatially continuous ODE represen- 
tation. We then incorporate a homogenized version of the 
individual induction motor dynamics to create a PDE repre- 
sentation of electrical feeder dynamics. 



A. Discrete Power Flow Model-Dist Flow Equations 

The flow of electric power in the quasi- static approxi- 
mation i s con trolled by the Kirchoff laws. The DistFlow 
equations 21 22 are these equations, written in terms of power 
flows and in a convenient form for the radial or tree-like distri- 
bution circuit with a discrete set of loads shown schematically 
in Fig. [4^, 



Pn+l — Pn 



u n+l 



-2(r n p n + x n (f> n ) 

- (r 2 +a; 2 ) P» + ^ 



(7) 
(8) 

(9) 



Here, n = 0, • • • , N — 1 enumerates the sequentially- 
connected buses of the circuit, and p n , (j) n are the real and 
reactive power flowing from bus n to n + 1. v n is the bus 
voltage, while P n and Q n are the overall consumption of real 
and reactive powers by the discrete load at bus n. The values 
of r n and x n are the resistance and reactance of the discrete 
line element connecting n and n + 1 buses. The voltage at 
the beginning of the line is nominally fixed by control equip- 
ment, and there can be no flow of real or reactive power out 
of the end of the circuit. These two observations provide the 
following boundary conditions: 



vo, Pn+i 0iv+i = 0. 



(10) 



Eqs. ( |7|8|9|10| ) combined with the given real and reactive con- 
sumption pattern, P n , Q n for n = 1, • • • , N, uniquely define 
profile of voltage, v n , and power flows, p n , <j> n9 along the cir- 
cuit. 
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where the conversion to continuous form consists of replacing 
X m , R m and P, Q, To, M by the respective densities x m , r m 
andp, g, to, and /i. The new boundary conditions are 

v(0) = i7 , = = 0. (17) 



^l™ ../tx Equations dll|17b form our PDE model of a distribution 



p(z),$Kz) feeder loaded with induction motors. We assume that the dis- 



tributions of all the density parameters along the circuit are 
known, and in this initial work, we assume these densities are 
constant. Note that evolution in the model occurs solely due 
FIG. 4: A feeder line modeled with a) a discrete set of electrical loads to temporal derivatives in Eq. ( [T4] ) representing mechanical 
and with b) a continuous distribution of loads. relaxation of the spatial distribution of motors. 



We note that the dynamical relaxation of distribution cir- 
cuit power flows p n and <p n will also occur on electrical time 
scales, i.e. much faster than the mechanical dynamics in 
Eq. ([5]). Therefore, the quasi-static pow er flow s described in 
the DistFlow formulatiorPKS m Eq S ( 7|8|9 ) is a sufficient 
starting point for the phenomena discussed in the manuscript. 



B. Continuous Power Flow Model 



IV. PHASE TRANSITIONS AND HYSTERESIS IN A 
FEEDER 

In this Section, we discuss the physical picture of phase 
transitions and hysteresis that emerges from analysis and sim- 
ulations of the 1+1 space-time PDE model of Eqs. ([TTJT7 ). In 
Section IV A we begin with a qualitative description of the 
hysteretic, phase transition-like behavior of the distribution 
feeder. Section lTVBl discusses numerical results in the frame- 



^ « . rra . . r ^ , ^. ™ . work of the qualitative arguments of Section IV A Then, in 

Following , the continuous form of the DistFlow equations c nrrT^i • * -a i — z • a- 



is 



d z p = -p-r- 



9 ' 

,2 , A2 



d z (j) = -q - x- 



(11) 
(12) 



Section [rVC| we regress again to provide general physics dis- 
cussion of the phase transition special features observed in the 
simulations. 



A. From Local (Motor) to Global (Feeder): Qualitative Picture 



where z is the coordinate along the distribution circuit, r, x 
are the per-unit-length resistance and reactance densities of 
the lines (assumed independent of z) and p(z) and q(z) are 
the local densities of real and reactive powers consumed by 
the density of the spatially continuous distribution of motors 16 
at the position z G [0; L] . The power flows p and <p are related 
to the voltage at the same position according to 16 



d z v 



rp + xcj) 



(13) 



C. PDE Model of Feeder Dynamics 

Instead of the standard voltage-independent (p, q) model of 
distributed loads, discussed in 16 , we consider the more com- 
plex dynamical loads described above. The load densities 
p(z) and q(z) are related to uj{z) and v(z) through the density 
versions of Eqs. ( 1|5|6 ) 



d 



p_ 



t 



sr m v* 



% + S Z X: 
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+ s z x, 



(14) 
(15) 
(16) 



From the qualitative description of distribution circuit and 
induction motor dynamics and FIDVR events given in Sec- 
tion [TJ we make an analogy between FIDVR and a first-order 
phase transition 26 where the z-dependent motor frequency 
lj(z) plays the role of the order parameter, i.e. u ~ uo in 
the "normal phase" and uj « in the "stalled phase". Re- 
casting Eq. ( [T4] ) in terms of an cj-dependent potential U, we 
find 



duo 
~dt 

where U — 



dU 
duo ' 



-I 

^0/ 



p(u';v) 



UJ 



duo' 



(18) 
(19) 



and p(uj;v) is defined by Eq. (15). The resulting effective 
potential U(uj;v) is shown graphically in Fig. [5] for different 
values of v. 

We identify three different regimes in Fig. |5j (a) high volt- 
age v > v+ (purple curve) where the normal state is the only 



stable solution, (b) intermediate voltage v c < v < 



(red 



curve) where the normal state and the stalled state coexist, i.e. 
motors can be in either of the two states, and (c) low voltage 
v < v~ (dark blue curve) where the stalled state is the only 
stable solution. The boundaries between these regions occur 
at the spinodal voltages v+ (light blue curve) where the state 
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FIG. 5: Local potential U(u) for different values of the voltage v and 
a = 0.1 demonstrating the continuous evolution of the normal and 
stalled states from only the normal state at high voltage (purple), co- 
existing normal and stalled states at intermediate voltages (red), and 
only the stalled state at low voltage (dark blue). The light blue curve 
is for the spinodal voltage - the boundary between the high and 
intermediate voltage cases. The green curve is for spinodal voltage 
v~ - the boundary between the low and intermediate voltage cases. 
The inset is a view of the region near uj — showing the emergence 
of a stalled uj « minimum of the potential for v < . The states 
labelled (a, 6, c, d) are the same as those in Figs.|l|3| 



(green curve) where 



jumps from c to a in Figs. [1|3] and at v t 
the state jumps from b to d in Figs. T]|T 



The high- voltage case (a) where the motors can only be in 
the normal state is the desired regime for electrical grid oper- 
ations. The low- voltage case (c) is undesirable and is a result 
of the electrical torque at low voltage not being able to over- 
come the mechanical torque with the subsequent decline of uj 
to very small values. Case (b) also allows for the motors to 
be in the normal state, but the existence of two minima in the 
intermediate voltage range v~ < v < v+ (see Fig. [5} is as- 
sociated with the overlapping high and low voltage states in 
Fig. [2] which can lead to local hysteretic behavior. For exam- 
ple, a small voltage perturbation can kick a motor from the 
normal state over the potential barrier (see red curve in Fig. [5} 
where Eq. ( [T8] ) shows that it subsequently relaxes to the stalled 
state. A simple reversal of the perturbation does not necessar- 
ily lead to the reverse transition. This entirely local hysteresis 
is important because it defines the possible states of motor 
operation, however, it is the long range interactions between 
the local motor behavior that creates phase transition-like be- 
havior with a phase boundary between separate normal and 
stalled phases. 

When a motor undergoes the local transition from a normal 
to a stalled state, its rotational frequency uj changes signifi- 
cantly. Motors distributed along the circuit do not interact via 
uj, however, changes in uj drive large changes in local reactive 
power density q (see Fig. [3]) which couples to all of the other 
motors via the power flows (p, <p) and voltages in Eqs. ( [TT} 
13 ). Crucially, if a perturbation causes a group of motors in 
the tail segment of the feeder (z ~ L) to enter the stalled 
states (uj « or s « 1), the increase in the local reactive load 



feeder, and Eq. ( [13] ) shows that v will be depressed at all z 
along the feeder. The voltage depression will be the largest at 
the tail (z ~ L), and if the depression is severe enough, the 
terminal voltage of the normal- state motors neighboring the 
stalled tail section will drop below v~ and the local potential 
U (uj) changes from the purple or light blue curves of Fig. [5] 
to the green or dark blue curves. Equation (18) shows that 
these motors relax into the stalled state, further increasing the 
power flows (j) everywhere along feeder. This phase transition 
front continues to propagate toward z = as the increases in 
local q drive increases in the non-local <j) which depress the 
non-local v. The voltage boundary condition at the head of 



the feeder (Eq. ( [17] )) may stop the the front before its reaches 
all the way back to z = 0, however, the establishment of a 
globally stalled or partially stalled phase is hysteretic because 
simply reversing the original spatially-local perturbation, i.e. 
returning the small section of motors in the tail to the running 
state, will not be able to overcome the globally-stalled state 
once it has become established. 



B. From Local to Global: Numerical Experiments 



Next, we perform numerical simulations of Eqs. (TTJT7]) to 
explore the qualitative dynamical description given above. We 
examine how one property of the feeder, i.e. its length L, 
affects the qualitative description given above by performing 
two identical simulations except L = 0.45 in Fig. [6] and is 
only slightly longer at L = 0.5 in Fig. [7] Although the change 
in length is relatively minor, the final states of the two feeders 
are radically different. In both simulations, we start with all 
the induction motors (the only type of load considered here) 
in the normal state, i.e. uj ~ ujo. At t = 0, a perturbation 
is applied where is abruptly lowered to 0.8 and held there 
long enough so that all of the motors stall. Subsequently, vq 
is restored to 1 .0 and the evolution of the motors and feeder 
variables are monitored. The forced evolution of emulates 
the behavior that would be driven by a nearby fault on the 
transmission system supplying the feeder and its substation. 

Figure [6^ shows the state of the L = 0.45 feeder immedi- 
ately after the fault is applied. The voltage (black curve) starts 
at v(0) = 0.8 and droops only slightly. Although v < v~ , the 
inertia of the motors (although small) keeps them temporarily 
at uj/ujo ~ 1 (blue curve) so that their local reactive loads q 
(pink curve) remain low as does the non-local reactive power 
flow (red curve). However, because v < v~ the motors no 
longer have an equilibrium state at uj/ujq ~ 1 (see dark blue 
curve in Fig. [5} and Eq. ( [T8] ) forces them to relax to the stalled 
state which is evident in Fig. [6]} where uj/ujq ~ (blue curve) 
all along the feeder. All of the induction motors have now 
made the transition to the upper branch of the reactive power 
curve in Fig. [3] near to state a which is reflected by the in- 
crease in q (pink curve) and the reactive power flow (j) (red 
curve). Subsequent to Fig. [6]), is restored to 1.0, and a left- 
to-right propagating phase front is formed (see blue curve in 
Fig. [6]:) where the normal and stalled phases are segregated 
with all the motors to the left of the front in the normal state 
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(a) (b) (c) 



FIG. 6: Sequence of snapshots of a simulation for a feeder with L = 0.45 < L c that undergoes a short voltage fault. In this case L is 
short enough so that the long-range interactions via power flows and voltage are insufficient to promote the local hysteresis at each motor into 
globally hysteretic behavior. Snapshot (a) shows the beginning of the sequence (immediately post fault) where the small inertia of the motors 
maintains cj/cjo ~ 1 (blue) so that the the reactive power load density q (pink) and reactive power flow <j) (red) are low and the voltage profile 
(black) relatively flat, although near 0.8. Snapshot (b) shows the feeder at a time after the reduction of electrical torque has decelerated all of 
motors to uj/ujq ~ (blue) with an associated rise in reactive load density q (pink) and reactive power flow (red) causing the voltage profile 
(black) to droop more than in (a). Snapshot (c) shows the situation after the fault is cleared with vo restored to 1.0. Although the reactive 
power flow 4> is still high, it is insufficient to force v < and the motors relax back to the normal state as the phase front propagates from 
left to right. The dynamical nature of this transition is evident from the real power load density p (yellow). The peak in p at the front is a 
result of the acceleration of the motors as they transition from the stalled to normal state (see Eq. (T4|). The motors shown accelerating in (c) 
eventually reach cj/cjo ~ 1, and the feeder relaxes to a globally-normal state. See Appendix [5] and Supplementary Information for respective 
movie (Movie 4: movie_recovery.pdf). 



while those to the right are in the stalled state. The motion 
of the finite-thickness phase front is evident from the local 
real power load density (yellow curve in Fig. [6]:). The peak 
in p (above the steady-state values observed far to the left or 
right of the front) is responsible for the acceleration of the mo- 
tors within the front, and as the motors in front accelerate to 
uj/ujo ~ 1, the front advances to the right. A left-propagating 
front would show a downward peak in p. After vo has been 
restored to 1.0, the reactive power q drawn by the stalled mo- 
tors along this feeder is insufficient to create a large enough 
(j) to lower v into the intermediate range v~ < v < v+, and 
the feeder completely recovers to its initial state, i.e. the re- 
versal of the initial perturbation restores the feeder state. Al- 
ternatively, the implicit long-range power flow interactions, 
expressed in a globally depressed voltage profile, are insuffi- 
cient to turn the local hysteresis into global hysteresis. The dy- 
namic version of this simulation is provided in the SI (movie 
4: movie_recovery.pdf). 

In Fig. [7] we show the final steady state of the exact same 
simulation as in Fig. [6] except that we have slightly increased 
the feeder length, L = 0.50. After restoration of vq = 1.0, 
phase front (blue curve) still forms and propagates into the 
feeder, however, it becomes stationary at z ~ 0.22. The lack 
of a peak in the real power load density (yellow) shows that 
there is no motor acceleration in the front implying that it is 
stationary. In this case, the reversal of the initial perturba- 
tion does not restore the original state and the feeder displays 
significant hysteresis. The added reactive power load density 
q between z = 0.45 and 0.50 interacts with the impedance 
over the entire length of the feeder to create conditions (i.e. 
v < v+) near z ~ 0.25 that enable the local motor hys- 
teresis to make the feeder globally hysteretic. The dynamic 



version of this simulation is provided in the SI (Movie 3: 
movie_hysteresis.pdf). We note that a large section of motors 
to the right of the stationary phase front have v~ < v < v+ 
and therefore also have a stable normal state (see Figs. [2] and 
[5]). Additional small perturbations could result in local transi- 
tions to the normal state and subsequent global recovery. 

In between L = 0.45 and 0.50 is a critical length L c where 
the hysteresis first appears. Feeders operating with L > L c 
can be called "dangerous" because they are operating nor- 
mally until a voltage fault occurs. Post fault, a fraction of the 
feeder does not recover to normal operation and this fraction 
of induction motors remains stalled, i.e. the feeder has just 
undergone a FIDVR event. Over a period of one to several 
minutes, the stalled motors will get disconnected by tripping 
of thermal protection systems, however, this uncontrolled re- 
covery may also lead to overvoltages that are equally trouble- 
some. 



C. Special Features of the Phase Transition 

Although the dynamical behavior described above is remi- 
niscent of many of phase transitions in physics, it also displays 
some very unique features, e.g. its lack of explicit spatial lo- 
cality and the instantaneous nature of the global voltage ad- 
justment. In the Ginzburg-Landau (GL) theory of first-order 
phase transitions, we also have a potential with two minima, 
however, the spatio-temporal dynamics of the phase transition 
are different. In contrast to the spatio-temporal dynamics of 



the current problem (as given in Eq. ([18])), GL dynamics are 
typically driven by a dispersion term, i.e. by adding a term 
such as d\u to the right hand side of Eq. ( 18 ) which effectively 
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FIG. 7: Voltage fault and subsequent hysteretic recovery for L — 
0.50 > L c . The sequence of events (fault, transient and initial re- 
covery) are identical to the one shown in Fig. [6] however the recov- 
ery is not complete. The additional loading from the motors between 
z — 0.45 and z — 0.50 generate long-range interactions that force 
v < v+ for z > 0.22. Although vo has been restored to 1.0, the mo- 
tors with v < do not recover to a normal state, even though many 
of them do have such a stable state. In this case, the long-range inter- 
actions are are sufficiently strong to promote the local hysteretic be- 
havior of each motor into globally hysteretic behavior. See Appendix 
[B]and Supplementary Information for respective movie (Movie 3: 
movie_hysteresis.pdf). 



phase-phase interface is subject to a boundary condition that 
relates its speed to the temperature at the interface. If the in- 
terface progresses, heat is released locally and it is then trans- 
ported via diffusion. Different versions of the problem show 
many interesting behaviors. In a semi-infinite domain, self- 
similar continually slowing fronts emerge. In finite domains, 
the fronts may become stationary. This behavior is similar to 
what we observe in the present problem - the analogy with the 
Stefan problem is the emergence of a global solution with an 
inhomogeneous order parameter profile and sharp interfacial 
boundary. However, the physics and interplay of the mecha- 
nisms that creates the behavior in the two problems are sig- 
nificantly different - thermal diffusion vs heat release in the 
Stefan problem, and local frequency transformations vs volt- 
age profile and rearrangement. Moreover, the latter point also 
emphasizes the difference - voltage profile plays the role of 
the heat injection but it does it globally along the feeder, also 
changing instantaneously, i.e. voltage rearrangement takes 
place with an infinite speed (electro-dynamic effects are in- 
stantaneous in our electro-mechanical model) while the heat 
injection condition is modified gracefully in the Stefan case, 
as the heat propagates along the domain with a finite speed 
controlled by thermal diffusion. 



V. SIMULATING AND EXPLAINING FAULT INDUCED 
DELAYED VOLTAGE RECOVERY (FIDVR) 



couples the order parameter uj at different spatial locations. 
The resulting GL dynamics is a phase transition front with a 
shape and speed defined by the local balance (i.e. within the 
front) of the added dispersion term against the existing non- 
linearity (rhs of Eq. ([T8])) and dynamics (lhs of Eq. dT8])). In 



contrast, the transition front dynamics of the present problem 
is driven by the globally-superimposed spatial inhomogene- 
ity of the voltage profile v(z). During short periods of time 
when there are not abrupt global changes, the voltage pro- 
file remains relatively frozen and the motors respond locally 
(Eqs. ( |18|19 )) to the mismatch between their current rotational 
state u and minimum U(u) as determined by the local voltage 
v(z). These adjustments occur most rapidly in the vicinity of 
the phase front where the state mismatch is the largest. The 
adjustment of v(z) to the evolving motor loads is instanta- 
neous (Eqs. ( pT]|13|15|17) ), but the evolution is actually tem- 
porally slow because the motor states are only adjusting in the 
small region of the phase transition front. 

Another point of comparison is the so-called Stefan prob- 
lem (se d 27 l 28 l and references therein) describing a phase tran- 
sition driven by heat released at the interface between the 
two phases. In its one dimensional formulation, the Stefan 
problem considers two sub-domains with their outer bound- 
aries (the boundaries away from the phase front) maintained 
at different conditions, e.g. one at a constant temperature flux 
and another at different temperature. Within each of the sub- 
domains temperature plays the role of the order parameter, 
and it obeys simple thermal diffusion (possibly with differ- 
ent diffusion coefficients in the two sub-domains). The sharp 



In the previous Sections, we built up an understanding of 
the different types of dynamical behavior of an induction 
motor-loaded distribution feeder and the different final states 
that may result. In this Section, we discuss in more detail 
the processes specific to FIDVR. Although we discussed the 
anatomy of a FIDVR event in the introduction, we repeat it 
here to motivate the following study of FIDVR dynamics. 
Prior to the transmission fault, the voltages vq at the head of all 
the distribution feeders extending from the substation served 
by the transmission line are 1.0, and all the feeders are in the 
normal phase, i.e. all the induction motors are in the normal 
state. During a transmission fault, the large fault currents in 
the transmission lines locally depress the transmission voltage 
which in turn depresses vq. vq remains depressed below 1.0 
until the transmission fault has cleared, i.e. automatic protec- 
tion circuit breakers have opened to de-energize the faulted 
line to extinguish the ionized air supporting the fault current 
and then reclosed to re-energize the transmission line. During 
the fault-clearing process, one transmission line (out of typ- 
ically two or more) serving the substation is briefly removed 
from service so which also tends to depress vq. It is during this 
period of vq depression that the induction motors on the dis- 
tribution feeders undergo "collapse dynamics" and may stall. 
If the transmission fault is cleared normally, the voltage at the 
the substation returns to near normal levels after the circuit 
breakers reclose. If there is significant motor stalling, v$ may 
not fully recover to 1.0 , but we approximate the post-fault 
voltage as vo = 1.0. The induction motor-loaded feeders then 
undergo "recovery dynamics" which evolve to a final steady 
state. 
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In Section [VA] we first discuss the collapse dynamics dur- 
ing a period of the fault when vo is depressed. This is followed 
by Section |VB| where we discuss the recovery dynamics (or 



partial recovery) after vq is restored. Finally, in Section V C 



we analyze the two periods in combination asking the ques- 
tion: under which conditions does a fault result in a FIDVR 
event? 



shape propagating with the speed and shape controlled by in- 
stantaneous voltage at the point of the front which defines the 
slope of the energy landscape U (cj) in Eq. ( 19 ) and Fig. [5] 



B. Post-Fault Dynamics-Recovery from Stalled (or Partially 
Stalled) States 



A. Dynamics During a Fault-Collapse to Stalled Sates 

If the voltage dr op dur ing the fault is large enough, i.e. if 

< v~ (see Figs. |l|2|5| ), all the motors will rapidly deceler- 
ate to uj/uo « and the entire feeder will eventually end up in 
the stalled state. Figure [8] (and movie l:movie_large _fault.pdf 
in SI) illustrates the details of this process which are consis- 
tent with our interpretation in terms of an "overdamped fall" 
down the the v < v~ energy landscape of Fig. [5] (dark blue 
curve) governed by Eq. ( [T8] ). We note that the long-range spa- 
tial coupling forces the motors farther down the feeder to col- 
lapse to uj/ujq ~ faster and reinforces the collapse after 
is starts. Specifically, the cumulative effect of the induction 
motors loads reduces the voltage at the remote locations of 
the feeder (black curve) resulting in a steeper local potential 
U (oj) (dark blue curve in Fig. [5]). Eq. ( [18] ) shows that these 
remote motors stall sooner, which is borne out in Fig. [5]: (blue 
curve). Additionally, as the remote motors decelerate, their 
local reactive load density q (pink) increases which increases 
(j) (red) further lowering v(z) (black) and increasing the slope 
of the local potential U(u). The result is an even faster col- 
lapse into a pure stalled phase. Moreover, we observe that the 
more severe the voltage drop or the lower the motor rotational 
inertia /i, the earlier these motors will be stalled. 

Next we consider a less severe fault where v~ < vq < 
1.0 (see Fig. [9] and the full movie version, Movie 2: 
movie_small_fault.pdf in the SI). Although vo > v~ , the im- 
mediate post-fault voltage drops along the feeder to a point zo 
where v(zq) < v~ . If the voltage profile was subsequently 
frozen in time, we would expect the motors with z > zq to 
behave very much as in Fig. [8] i.e. a decleration to uo/uoq ~ 0, 
although somewhat slower in the previous case because the 
local v(z) is slightly higher. In Fig. J9J), we do observe this 
initial behavior, however, the local reactive power loading q 
(pink) again increases as the motors stall increasing </> (red) 
and lowering the local v(z). Via this spatial coupling, the 
feeder power flows reinforce the collapsing wave allowing it 
to propagate to z < zo . 

Although the reinforcement process can be significant, it 
cannot overcome the the boundary condition at z = that 
maintains vo > v~ . As shown in Fig. [9jl, the phase front fi- 
nally stops at location < z < zo where the local voltage 
finally stabilizes at v(z) = v~ . As this occurs, local motor 
accelerations cease and the phase front in oj, p, and q steepen 
creating a sharp demarcation between the two pure phases of 
normal and stalled states. One common and possibly univer- 
sal feature in the two scenarios discussed above is emergence 
of a "collapse" transient which can be interpreted as a quasi- 
stationary phase transition front of a slowly evolving "soliton" 



The severity of the vo depression during the fault deter- 
mines the final state that will be reached if the collapse dy- 
namics are allowed to evolve to steady state. The transmis- 
sion fault may be cleared (vq returns to 1 .0) before this steady 
state is reached, and we explore the dependence of the recov- 



ery dynamics on the fault time duration in Section V C For 
simplicity of the present discussion, here we assume we are 
starting from a post-collapse steady state. Whether the feeder 
starts in a partially or fully-stalled post-collapse steady state, 
the perturbation of v$ returning abruptly to 1.0 has the chance 
of restoring the feeder to the fully normal phase. 

We have performed a range of simulations starting from 
both partially and fully-stalled steady states resulting from 
simulation of a depressed of vq. We then restore vq = 1.0 and 
monitor the recovery dynamics. One such case for a fully- 
stalled initial condition that fully recovers is shown in Fig.[T0| 
When restoring vq = 1.0 is sufficient to drive a full recovery 
(as in Fig. [TO} , we observe a rather rich dynamics which can 
be split, roughly, into the following three stages: 



Figures [TOp— ^b: The abrupt change of vq instanta- 
neously creates a smooth voltage profile dependent 
upon the instantaneous adjustment of real p and reac- 
tive q load densities to the higher voltage, but at their 
motor's stalled rotation speed. Motors with v(z) > v+ 
start to accelerate, and the phase transition front is set 
up in the vicinity of the feeder head. 



• Figures [T0p—^e: The recovery front matures and sharp- 
ens while propagating into the feeder. The narrow re- 
covery front is always located near to v(z) = v+ where 
the fast dynamics of motor acceleration and state tran- 
sition occur. These fast transitions are accompanied by 
jumps in p and q which act in a non-local manner to 
push the location of v(z) = v+ to larger z thus driving 
the transition front farther into the feeder. The propa- 
gation, although driven by the fast dynamics of motor 
state transitions, comprises a slow dynamics because 
the state transitions occur in a phase front that is narrow 
compared to the feeder length. Also part of the phase 
front is a peak in p above the steady- state values on ei- 
ther side of the front. This peak is required to supply the 
power to accelerate the motors from a stalled to a nor- 
mal state. If the feeder is long enough, the phase front 
appears to reach a time-invariant shape that propagates 
much like a soliton. 



• Figures [T0j—^h: The slow dynamics of the recovery ap- 
proaches to within about one or two phase front widths 
of the end of feeder. With fewer motors to acceler- 
ate, the voltage adjusts faster and the front propagation 
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FIG. 8: Sequence of snapshots illustrating the collapse dynamics during a large voltage fault with vo < v~ and L — 0.6 > L c . Snapshot 
(a) shows the situation just after the application of the fault - although the v(z) < v~, the motors continue to rotate with uj/ujq ~ 1 because 
of their inertia. Snapshot (b) corresponds to a short time after the application of the fault when the motors are just starting to decelerate. The 
more remote motors experience smaller v and steeper U(uS). Their faster deceleration results in smaller uj /ujq at these remote locations. The 
collapse is reinforced by an increase in q (pink) as the motors reach lower u. Snapshot (c) is taken later in the process: the end of the line is 
completely stalled, the wave of deceleration starts to propagate backwards, from the tail to the head. Snapshot (d) shows the final phase: all 
the motors are stalled. Note a significant increase in the reactive power drawn by the stalled feeder. See Movie 1 : movieJarge_fault.pdf in SI 
for a dynamic version of this process. 



L=0.6, time=0.46 



L=0.6,time=1.48 




0.2 0.4 0.6 

position along the feeder z 

(a) 



0.2 0.4 0.6 

position along the feeder z 

(b) 



0.2 0.4 o.e 

position along the feeder z 

(c) 




0.2 0.4 0.6 

position along the feeder z 

(d) 



FIG. 9: Collapse dynamics caused by a small voltage fault with vo > v~ and L = 0.6 > L c . The collapse proceeds slower than in Fig. [5] 
because the higher voltages result in shallower U(uj). Snapshot (a) corresponds to a short time after the fault as the motors have just begun to 
decelerate, however, the rise in q (pink) is already reinforcing the collapse. Snapshot (b) corresponds to later during the fault where the end of 
the line is stalled and the deceleration front is starting to feel the influence of the boundary condition at z = 0. Snapshot (c) the deceleration 
front has become nearly stationary as the long-range interactions can no longer overcome the boundary condition at z = 0. Snapshot (d) show 
the feeder in the stabilized, partially- stalled steady state. See Movie 2: movie_small_fault.pdf in SI for a dynamic version of this process. 



speed increases until is has consumed the entire feeder 
and the feeder reaches a uniform normal phase. See Fig. 
[TOi-h for illustration. 

See Movie 4: movie_recovery.pdf in SI for the full movie of 
the recovery phenomenon. 

Perhaps the most remarkable feature of the recovery pro- 
cess is formation, at the intermediate stage, of a "soliton" - 
a quasi-stationary shape moving with roughly constant speed 
from the head to the tail. The emergence of the "soliton" is 
due to the long-range interaction of the fast but local changes 
in p in q in the phase front to the slower but global changes in 
voltage v(z). 

As we will see in the following Section, the feeder does not 
always fully recover simply because vo is restored to 1.0. In 
these cases, we have explored the effects of temporarily (or 
permanently) raising vo above 1.0 and found this approach to 
be quite effective in restoring feeders that would have other- 
wise remained stalled. However, such intelligent control ac- 
tion requires reliable detection of the entry into a FIDVR event 
and fast switching and/or device control to increase vq to a 



sufficient level. We postpone full analysis of such a control to 
future work. 



C. Dynamic Transition: Will a Feeder Enter a FIDVR State? 

Our simulation results suggest that the normal phase, the 
fully stalled phase, and any of the partially stalled phases (with 
the feeder split in the normal head and stalled tail) can be the 
final stationary and stable point of a dynamical evolution. In 
this Section, we explore the properties of the final state as the 
properties of the feeder (length L and motor inertia /i) and 
the fault (magnitude of voltage depression Av and duration 
) are varied. Our goal is to develop an initial under- 
standing of the "non-equilibrium phase diagram" that controls 
whether or not the feeder recovers to a fully normal phase. In 
the test discussed below, we consider a feeder with L > L c 
because those with L < L c are known to always recover to a 
fully normal phase no matter the size or duration of the pertur- 
bation. In our study of the phase diagram, we dissect the four- 
parameter space, (L, /i, Av, T pertu ) in six different ways: 




FIG. 10: Eight snapshots of the evolution of a feeder from a steady-state, fully-stalled phase to the fully-normal phase following the restoration 
of vo = 1.0. (a) The final fully-stalled steady state reached after simulating the feeder for avo = 0.8 fault voltage depression, (b) vo is restored 
to 1.0 and the motors at the head of the feeder start to accelerate, (c) A recovery front is built which starts to propagate into the feeder, (d) The 
front continues to propagate at roughly constant speed and also showing a universal "soliton"-like shape, (e) The front continues to propagate 
with roughly the same speed and a universal "soliton"-hke shape, (f) Completion of the recovery process where interaction with the end of the 
feeder accelerates the recovery front, (g) The recovery front still propagates while the shape of the front starts to change because fo interactions 
with the end of the feeder, (h) End of the recovery process: the entire feeder is back to the normal state. See Movie 4: movie_recovery.pdf in 
SI for dynamical illustration. 



1. fixing ii and T pertu and exploring the (L, Av) subspace, 
see Fig.fTTj 

2. fixing fi and Av and exploring the (L, T pertu ) subspace, 



see Fig. 12 



3. fixing Av and T pertu and exploring the (L, ji) subspace, 



see Fig. 13 



4. fixing L and fi and exploring the (T pertui Av) subspace, 



see Fig. 14 



5. fixing the values of L > L c and (sufficiently large) A v 
and exploring the (/i, T pertu ) subspace, see Fig. 
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6. fixing L and T pertu and exploring the (/i, Av) subspace, 



see Fig. 16 



Starting as usual with the feeder in a normal state with vo = 1, 
we apply a vo depression of magnitude Av and duration 
T pe rtu with subsequent recovery to vo = 1 and then integrate 
the dynamics until the feeder reaches a steady state. Unless 
all of the motors on the feeder recover to a normal state, the 
feeder final state is classified as partially stalled. In all of the 
subsequent Figures, the filled red circles indicate this bound- 
ary between a fully normal phase feeder and a partially- stalled 
feeder. 

From our studies of the subspaces and the following Fig- 
ures, we can make the following conclusions: 



• From Figs. [TT|[T3j it is apparent that the feeder becomes 
less resilient to perturbations as it grows in length and 
that there is an upper feeder length L* « 0.63 that 
requires an infinitesimal perturbation to force it into 
a partially-stalled phase. In fact, if the feeder is too 
long (i.e. L > L*), the normal phase is no longer 
stable (v(L) < v~) and the feeder always has a par- 
tially stalled phase. Therefore, for the remainder of this 
study, feeder with L c < L < L* are of interest because 
feeders with L outside this range are either robust to 
all perturbations or always unstable to a partially- stalled 
phase. 



• From Fig. [14] it is apparent that there is value of 
Av ^0.05 such that the feeder will never stall no mat- 
ter how long the perturbation is applied. At Av = 0.05, 
vo = 0.95 and the static voltage drop at these reduced 
voltages would make v(L) « 0.83 = v~ . From this in- 
terpretation, this lower bound on the Av for extremely 
long T pertu can be computed from the static power flow 
equations by looking for the vq that forces v(L) = v~ . 

• The two previous conditions can be computed from 
static considerations. The interesting dynamics be- 
hind the transition to a partially-stalled state is then ex- 
pressed in Figs. [15] and [16] which can be understood 
by considering the decline of the motor rotational fre- 
quency that occurs at the far end of the feeder dur- 
ing the time T pertu of the fault. We can crudely ap- 
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FIG. 1 1 : Phase diagram of the dynamic transition from the normal 
phase to a partially stalled phase. All the points above the curve lead 
to a partially stalled phase, all the points under the curve result in a 
normal state. For L < L c ~ 0.46, there is no hysteresis and the 
system always ends in a fully running phase. 



FIG. 12: Phase diagram of the dynamic transition from the normal 
phase to a partially stalled phase. All the points above the curve lead 
to a partially stalled phase, all the points under the curve result in a 
normal phase. For L < L c ~ 0.46, there is no hysteresis and the 
system always ends in a fully running phase. 



proximate this decline by the product of fault dura- 
tion and the rate of frequency decline immediately after 
the application of the fault, i.e. A(uj / uj )\ t=T ertu w 



— (2P/vqLJq) (T pertu Av/[i) where we have used Eq. 14 
and a linear expansion of Eq.[l5] It is reasonable to ex 
pect that the boundary between a feeder that fully recov- 
ers and one that is partially stalled would be expressed 
by A(u/uo) \t=T pertu ~ constant. If such a relationship 
is found to hold, it would imply (T pertu Av / fi) « con- 
stant, which within scope of our parametric study, is in 
rough agreement with the results in Figs. [15] and [TSJ es- 
pecially if we account for the minimum required value 
of Av discussed immediately above. 

The approximate analysis above provides a good qualitative 
understanding of the fault and feeder parameters that lead to 
a feeder with a partially-stalled phase. However, more rigor- 
ous analytical analysis of Eqs. [TT][T7] is required to put these 
conclusions on firm footing. 

Beyond just determining the boundary between a fully- 
normal and partially-stalled feeder final states, the character- 
istics of the voltage fault also determine the number of mo- 
tors that will be stalled after the fault is cleared. The more 
severe the fault (in duration T pertu or in amplitude Av), the 
more motors will get stalled, up to a point. If the fault is se- 
vere enough (Av > Av* or T pertu > T* ertu ), the number of 
stalled motors will not increase any more. For example, the 



maximum number of stalled motors is reached in Fig. 17 i for 
a given set of parameters of the feeder line, and an even more 
severe fault will end up in this same state after it is cleared. 

On the other hand, the less severe the fault the less the 
number of stalled motors, and the system can reach a con- 
tinuous set of stable partially-stalled phases between the max- 
imum number of stalled motors and none of the motors stalled. 
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Fig. 17 shows examples of the different phases the system can 
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FIG. 13: Phase diagram of the dynamic transition from the normal 
phase to a partially stalled phase. All the points above the curve lead 
to a partially stalled phase, all the points under the curve result in a 
normal phase. For L < L c ~ 0.46, there is no hysteresis and the 
system always ends in a fully running phase. 



reach, and the file movie_different_states.pdf in SI illustrates 
the dynamics leading to these final states. 



VI. DISCUSSIONS & PATH FORWARD 

In this manuscript, we have modeled, simulated, analyzed 
and explained how the electro-mechanical dynamics of in- 
duction motor-loaded distribution feeders can lead to Fault- 
Induced Delayed Voltage Recovery (FIDVR), an effect ob- 
served recently in power distribution feeders. We approached 
these dynamics and the FIDVR events from the stand-point 
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FIG. 14: Phase diagram of the dynamic transition from the normal 
phase to a partially stalled phase. All the points above the curve lead 
to a partially stalled phase, all the points under the curve result in a 
normal phase. 



FIG. 16: Phase diagram of the dynamic transition from the normal 
state to a partially stalled phase. All the points above the curve lead 
to a partially stalled phase, all the points under the curve result in a 
normal phase. 
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FIG. 15: Phase diagram of the dynamic transition from the normal 
state to a partially stalled state. All points under the curve lead to a 
partially stalled regime, all points above the curve result in the nor- 
mal final state. 



of physics - explaining and interpreting them as an instance 
of a broader class of nonlinear electro-mechanical phenomena 
in power distribution systems. The main ideas and questions 
discussed in this work can be summarized as follows: 

• The 1+1 (space+time) continuous model, introduced 
irP^ and developed here, offers a computationally effi- 
cient framework for analysis of nonlinear and dynami- 
cal phenomena in distribution feeders, and crucially, the 
PDE model enables analogies with other known spatio- 
temporal dynamical systems and solutions to help build 
intuition about the dynamical behavior of the electrical 
grid's electro-mechanical dynamics. 

• By coupling a spatially local model of an individual 



motor - that describes its nonlinear, bi-stable and hys- 
teretic switching between normal and the stalled states - 
to a continuum version of the electrical power flow 
equations, we are able to explain the coexistence of a 
spatially-extended normal phase with a stalled phase 
(motors at the head portion of the feeder are in the nor- 
mal state, while motors in the tail portion of the feeder 
are stalled). 

• The emergence of the multiple spatially-extended states 
is interpreted in terms of a first-order phase transition 
where the distribution of the motor rotational frequency 
along the feeder is the order parameter. The voltage 
distribution along the feeder plays the role of an ex- 
ternal field (degree of freedom) that modifies an ef- 
fective energy potential for motor frequency. Different 
(normal, stalled or partially stalled) spatially-extended 
phases are stabilized and reach a steady state by achiev- 
ing a global electro-mechanical balance of their voltage 
and frequency distributions along the feeder. 

• Sufficiently strong perturbations, e.g. sudden drops 
or rises in voltage at the head of the feeder, lead 
to transients in the form of propagating phase fronts 
that separate normal-state motors at the head of the 
feeder from stalled-state motors in the tail of the feeder. 
We analyzed and classified the different types of tran- 
sients and resulting steady-state phase distributions that 
emerge after the transients settle for different perturba- 
tion strengths and lengths of the feeder. 

• We also experimented with the dynamics of mechanical 
frequency and voltage phase distributions under more 
realistic, but more complex, two-step perturbations - a 
voltage drop at the head of the feeder followed shortly 
by restoration back to its nominal value - that are ex- 
pected to simulate the real world perturbations that re- 
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(a)T pertu = 0.221 and Av = 0.15 (b)T pert u = 0.23 and Av = 0.15 (c)T pertu = 0.3 and Av = 0.15 (d)T pertu = 0.5 and Av = 0.15 



FIG. 17: Different final states can be achieved by manipulating the duration of the perturbation (T per tu), while its amplitude Av remains 
constant. From left to right, the fault is longer and longer and the number of stalled motors increases. See file movie_different_states.pdf in SI 
for five movies illustrating the different phases that can be reached and the dynamics that lead to the phases. 



suit in FIDVR events. The dynamics and emerging 
steady states were explored for different voltage pertur- 
bations (depth and duration of the voltage depression) 
and feeder characteristics (feeder length and inertia of 
the connected induction motors). 

Major conclusions drawn from our numerical experiments and 
analysis are 

• Hysteresis. The system is strongly hysteretic: reversing 
a perturbation does not lead to a simple reversal of the 
dynamical trajectory. 

• Recovery Conditions. When the feeder is short enough 
or the voltage perturbation is weak enough (small 
enough in amplitude or short enough duration), the 
feeder recovers to a fully normal phase following a volt- 
age perturbation thus avoiding a FIDVR event. Longer 
feeders or stronger voltage perturbations lead to incom- 
plete recovery and, by modifying the three parameters 
beyond the recovery threshold, one explores a continu- 
ous family of different partially- stalled phases. 

• Self-Similar Transients. When a feeder is sufficiently 
long, recovery transients appear to show universal 
soliton-like phase fronts with normal phase propagating 
into the stalled phase with an (approximately) constant 
speed and time-invariant shape. 

This manuscript opens up a new line of research into 
physics-based analysis of transients and phase transitions in 
distribution feeders. We plan to continue this work focusing 
on the following generalizations and extensions: 

• We modelled distribution feeders as consisting of iden- 
tical induction motor loads distributed uniformly along 
the feeder. In reality, the motors may be different and 
distributed non-uniformly along the feeder, their distri- 
bution and parameters may fluctuate in time, and they 
are present with many other different types of loads. We 
will extend our purely deterministic analysis to a proba- 
bilistic framework to describe the effects of these forms 
of disorder and noise to resolve questions such as: what 
is the probability that the feeder with a given level of 
disorder will recover after a perturbation of given am- 
plitude and duration and not enter a FIDVR state? 



• As recently shown irP^, controls associated with dis- 
tributed generation, e.g. inverters coupled to distributed 
photovoltaic (PV) generation, can also be incorpo- 
rated into the spatially continuous modelling frame- 
work (ODE framework). The static model of 29 sug- 
gests that a feeder with a sufficiently large penetra- 
tion of the PV generation may show a rather compli- 
cated bifurcation diagram with the emergence of mul- 
tiple low- voltage solutions that are related to the low- 
voltage, stalled solutions discussed in this manuscript. 
Much like the effects discussed here, these inverter- 
driven low- voltage states are not the result of the be- 
havior of a single inverter, but rather result from the 
collective action of many inverters and their and inter- 
action with the nonlinearity and nonlocal behavior of 
the power flow equations. To address these dynami- 
cal problems in a more comprehensive manner, we will 
extend dynamical description of this manuscript to the 
case of distributed generation, and more generally, to 
the case of a feeder containing a mixed portfolio of dis- 
tributed generation and different types of loads, includ- 
ing induction motors. 

• In our first physics-based paper on the subject of 
electro-mechanical dynamical transients in distribution 
feeders, we relied mainly on numerical experiments. 
However, the problem formulation may allow asymp- 
totic theoretical analysis, in particular: accurate res- 
olution of the phase transition boundaries, determina- 
tion of the bifurcation (spinodal) points, propagation 
of soliton-like phase fronts, and analysis of the tails of 
the distribution functions that account for the aforemen- 
tioned effects of noise and disorder. 

• We have focused primarily on describing electro- 
mechanical dynamics and the perturbations and tran- 
sients that lead to FIDVR events. Armed with the com- 
prehensive understanding gained in the process, we are 
ready to attack the larger question of distribution feeder 
voltage control. Specifically, what is the least control 
effort needed to avoid a FIDVR event following a given 
type of fault? 
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Appendix A: Methods 



1. Analysis of Static ODE 



To solve numerically Eqs. ( |11|12|13|14|15|16| ) in the sta- 
tionary (time-independent) case, we proceed as follows: 

• using the mechanical torque t = to (jj^j (i- e t £ 

[0, to])» as a scanning parameter, we "map" the values 
of p, q, uj and v through equations ( |14|15|16| )); 



• for all these "maps", we choose only one part of the 
hysteresis curve (either the one when we start from low 
voltage, or the one from high voltage); 



• we then solve equations ( |11|12|13| ) using the "maps" of 
the other variables and the matlab bvp4c solver. 

Note that some of the stationary solutions discovered in the 
results of dynamic exploration were missed by the static anal- 
ysis detailed above. 



2. Dynamic Simulations 



To solve numerically Eqs. ( |11|12|13|14|15|16 ), we employ 
space-time discretization and use an explicit finite differences 
scheme. Our dynamical simulations are split in two steps. 

• During the first sub- step we take p and q fixed (outputs 
of the previous time step) and solve the feeder-global 
static Eqs. ( 11|12|13 ) for v, (j) and p under conditions of 
the fixed voltage at the head of the line and zero fluxes 



at the end of the line (Eq. ([17])). This sub-step instanta- 
neously imposes a spatially smooth and globally corre- 
lated voltage profile. 

• Once the global, p, </> variables are fully updated, we 
start updating on the second sub- step the local variables, 
0J,p,q. For oj, we use its explicit time dependence: 



(Al) 



where A t = T/N t is the time step (N t is the number 
of steps, and T the time of the simulation). Once uo is 
updated, we easily get the updated values of p and q 
using Eqs. ( 15|16| ). 




Appendix B: Movies 

a. Movie 1 : "movie Jarge -fault.pdf " 

This movie (see^) shows the dynamical transient follow- 
ing a large drop in vq to 0.8, i.e. vq < v~ . At t < 0, the 
line is in a stable stationary state where all the motors in the 
normal state. At t = + (the first picture in the movie), vo 
suddenly drops from 1 to 0.8, and the voltage profile along 



the line (black curve) immediately responds. In response to 
the lower v(z),the motors' rotational frequency (blue) decel- 
erates everywhere along the line at a rate dependant on the 
local value of the voltage - motors closer to the head of line 
where the voltage is larger decelerate at a lower rate. Motors 
at the far end of the line decelerate faster and become stalled 
first. The stalling of these remote motors results in an increase 
in their local reactive load q (pink) and the overall reactive 
power flow <fi (red) which reinforces the reduction inv(z) cre- 
ating a normal-to-stalled transition front that propagates from 
the tail to the head of the feeder. Eventually, the entire feeder 
becomes fully stalled. 



b. Movie 2: "moviejsmalLfault.pdf" 

This movie (see^) shows the dynamical transient following 
a small drop in to 0.95, i.e. vq > v~ . The response of the 
voltage and the mechanical frequency of the motors in the tail 
portion of the line is similar to that in Movie 1 - the voltage 
near the feeder tail is too low and they decelerate with those 
nearer the tail getting to the stalled u/uq ~ state first. The 
transition to the stalled state expands as it is again reinforced 
by the increase in reactive loading q (pink) and reactive power 
flow cj) (red). However, the phase front slows down, sharp- 
ens, and eventually stops propagating as it begins to feel the 
boundary condition at z = that does not allow the motors 
near the head of the feeder to stall. The result is a half- stalled, 
mixed phase distribution with motors in the tail part of the 
line stalled, while motors near the feeder's head in the normal 
state. The slowing down of the front is related to the effect 
of critical slow down typical of first-order phase transitions as 
they approach a spinodal point, i.e. as v(z) approaches v~. 



c. Movie 3: "movieJiysteresis.pdf" ' 

This movie (see 32 ) shows hysteretic behavior of a feeder 
with L > L c , i.e. the application of an opposite perturba- 
tion does not lead to a reversal of the dynamical trajectory 
back to the original state. From t = to 0.9, the movie is 
the same as Movie 1 "movie_large_fault.pdf", i.e. vq is re- 
duced to 0.8 at t = and the feeder and shows a complete 
collapse into the fully stalled phase. At t = 0.9, the per- 
turbation is reversed as is restored to 1.0. The increased 
voltage causes the motors at the head of the feeder accelerate 
and return to the normal state and a propagating phase front is 
formed (more details about this recovery front can be found in 
the description of Movie 4 "movie_recovery.pdf"). The front 
advances towards the tail until it reaches the point where the 
voltage has fallen to to v+ (voltage above which a stalled mo- 
tors accelerates). At this point, the phase front can proceed 
no further and we end up with a partially stalled feeder. The 
voltage is unable to increase further because of the non-local 
reinforcement of the voltage drop by the stalled motors be- 
yond the point where v = v+. Similarly to what was seen 
in Movie 2 (and what is generally the case when the system 
approaches a spinodal point) the dynamics gets slower as the 
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moving front approaches the critical voltage i> + . 

d. Movie 4: " movie -recovery. pdf" 

This movie (see^) is identical to Movie 3, except that the 
length of the feeder line now is significantly shorter with 
L < L c . In this case, there is no hysteretic behavior after 
a voltage fault is cleared — the feeder recovers completely 
because the shorter L has diminished the non-local effects of 
the motors near the tail of the feeder. The recovery process 
take a long time allowing us to distinguish different phases of 
the process. First, at 1 < t < 1.3 the beginning of the line 
accelerates quickly; then, at 1.3 < t < 9.5, one identifies 
emergence of the recovery front, which propagates down the 
feeder with a nearly stationary soliton-like shape. The recov- 
ery front starts to change shape as it approaches and interacts 
with the boundary condition at the feeder tail. The absence of 
motors to accelerate beyond z = L leads to fast recovery of 
the normal states in the vicinity of the tail for t > 9.5. 

e. Movie 5: "movie -different states. pdf' 

There are in fact five movies (see^) glued together in one 
file with a different movie on each page of the .pdf . Each 
of the movies shows the dynamics of the feeder line during a 
voltage fault and after the fault is cleared. Characteristics of 



the feeder and faults are identical in all the movies, except we 
increase T pertu from segment to segment. We observe that 
the number of motors stalled in steady state increases with 
in T pertu showing that different partially- stalled states can be 
reached depending on the duration of the fault. However, we 
also notice existence of an upper limit to the number of motors 
that can be stalled, i.e. no matter how severe the fault, half of 
the line will always recover. The two last movies in the file il- 
lustrate this phenomenon, i.e. even though the fourth segment 
has T pertu =0.5 and the fifth segment has T pertu = 1, they 
both result in identical steady states with the same number of 
stalled motors. 
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